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ABSTRACT 



Deep water directional wave spectra, measured by a NDBC 
3 -meter buoy off Monterey Bay, are transformed to shallow 
water using the linear refraction model by Dobson (1967). 
The transformed directional spectra are compared with 
measured spectra using pressure gauge arrays in shallow 
water at Marina and Santa Cruz. Refraction (Kr) , shoaling 
(Ks) and Jacobian (J) transfer functions are computed. 

The modeling results and measured wave data at Marina 
are in good agreement for easterly waves at high frequen- 
cies. On the other hand, the linear refraction model gives 
highly variable values of wave solutions at Santa Cruz. It 
is believed that the unrealistic wave results at Santa Cruz 
are associated with the complex bathymetry of Monterey Bay. 
The present work examines the accuracy and limitations of 
linear refraction model by field observations. 
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I. 



INTRODUCTION 



Wave climatology plays an important role in coastal 
dynamics as most energy is put into coastal region by wind- 
generated waves. Initially, ocean wave analysis was limited 
to a single wave train defined by a wave height and a period 
in a one-dimensional spectrum, but in the last couple of 
decades the emphasis has been placed on two-dimensional 
spectra . 

The shoaling waves undergo refraction, reflection and 
diffraction mostly influenced by the wave approach angle and 
the bathymetry of the area. A spectral approach requires 
the superposition of a large number of progressive sinusoi- 
dal component waves. Therefore, waves of small amplitude 
(linear wave) are assumed. This assumption neglects 
reflection, diffraction and postulates long-crested waves. 

The area of Monterey Bay has been selected for this 
comparative study of ocean waves. This area presents an 
interesting problem for predicting directional wave spectra 
because of its variable shoreline orientation relative to 
prevalent eastward waves and the presence of the Monterey 
submarine canyon bisecting the bay. The Monterey submarine 
canyon is the largest in the western hemisphere and causes 
substantial perturbations on the nearshore wave field. 
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Historically, Monterey Bay was one of the first, and 
presumably the most intensively studied, locations for wave 
refraction. Wave refraction diagrams of Monterey Bay were 
first constructed in 1948 by Johnson, O'Brien and Isaac 
(Wiegel , 1964). The first intensive refraction studies were 
performed for amphibious landing exercises at Fort Ord 
(1954) followed by studies by the Corps of Engineers (1958) 
for small craft harbor at Santa Cruz. Wave refraction and 
littoral drift calculations have been performed by Dorman 
(1969) and Arnal et al . (1973). All these studies were 
performed using hand-drawn refraction diagrams. More 
recently, numerical refraction calculations have been made 
for a number of specific marine construction sites (e.g., 
Thornton (1978, 1980)) and erosion studies (Oradiwe, 1986; 
McGee, 1986) . 

The objective of this study is to compare ocean wave 
measurements with the linear refraction model. The deep 
water two-dimensional spectrum is transformed to shallow 
water spectrum using transfer function computed from a 
linear wave refraction model by Dobson (1967). The deep 
water spectrum, S 0 (fo/ a o) i- s determined from data acquired 
from the National Data Buoy Center (NDBC) using the Fourier 
analysis technique as proposed by Longuet-Higgins (1963) and 
Cartwright (1963). The transformed spectrum S(f,a), is 
compared with measured shallow wave spectra for which the 
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data had been acquired from pressure gauge arrays located at 
Marina and Santa Cruz. 

The background theory is described in Chapter II. The 
data acquisition and equipment are mentioned in Chapter III. 
Chapters IV and V contain discussions and the conclusions. 
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II . THEORETICAL BACKGROUND 



A. WAVE REFRACTION 

1 . General 

The refraction of ocean surface waves on the 
continental shelf is necessary for forecasting ocean waves 
in shallow water. Practical techniques for calculating 
refraction diagrams were developed by Johnson (1948) and 
Munk and Arthur (1952). The graphical technique developed 
by these investigators is called the Orthogonal Technique. 
With the advance of digital computers, the numerical methods 
were introduced by Griswold (1963), Wilson (1966) and Dobson 
(1967). The numerical calculations of wave orthogonals (or 
rays) can be quickly and accurately plotted for engineering 
purposes. The study of wave transformation that includes 
refraction over complex underwater topography and in large 
areas is more commonly accomplished through numerical 
models. (Hayes, 1980) 

2 . Refraction Theory 

In this study, linear wave theory is applied to 
shoaling and refracting waves, which neglects reflection and 
diffraction assuming energy propagates only in the direction 
in which the waves travel. It also neglects amplitude 
dispersion. The wave refraction is based on generalization 
of Snell's law for phase speed variation and conservation of 
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energy. The phase speed of the wave depends upon wave 
period (or frequency) and the depth of water at particular 
locations. Wave energy is conserved between wave orthog- 
onals, where orthogonal is perpendicular to the wave crest. 
The orthogonal indicates the direction in which shoaling 
wave progresses. Based on linear wave theory, it is assumed 
that the wave period will remain constant but the 
wavelength, the velocity and the height, will vary as the 
wave shoals toward the beach. 

In general, the phase speed of a small-amplitude 
linear water wave is expressed as 

C 2 = £ tanh(kd) (2.1) 

where g is acceleration of gravity, k the wave number and d 
is the depth of water at a particular location. Consequent- 
ly, in shallow water the phase speed of the wave varies with 
depth only and Snell's law becomes applicable in the case of 
straight and parallel contours as 



sin g 
sin ciq 




( 2 . 2 ) 



where a is an angle between the wave crest and local depth 
contours and C is phase speed. The subscript zero refers to 
deep water conditions. 
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The energy of a wave when it reaches the beach is a 
function of the wave height, and the wave height at the 
shore is related to the wave height in deep water. The 
variation of wave height within two wave rays is described 
(Mei, 1980) as 



H _ H r Cg 0i 1/2 r Vl/2 

H " H 0 [ Ti' ] [ b J 



(2.3) 



or 



H = H 0 K S K r 



(2.4) 



where 



K S = 

Kr = 



r Cg 0 -I 1/2 

1 Cg J 



_ 0 ! 1/2 
b J 



Shoaling Coefficient 



Refraction Coefficient 



(2.5) 



( 2 . 6 ) 



where b is the distance between two rays and Cg is group 
velocity and may also be represented in the form of 



Crr = 



= - [1 
2 1 



+ 



2 kd 



sinh 2kd 



Using the relationship C = C 0 tanh(kd) (Kinsman, 1965), 
Equation (2.5) may be written as 
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(2 . 5a) 



Ks ‘ V 



tl + 



2 kd 



sinh 2kd 



- 1/2 



From Equation (2.2), the change in angle determines the 
change in wave crest or b value. Thus the value of Kr may 
be computed from 



b 



0 



COS C*Q 



b 

cos a 



or 



% 



b 0 1/2 _ r COS a 0 1/2 
1 b J 1 cos a J 



(2 . 6a) 



where 



. -l r ^ sin Q 0 1 
a = sm [C — p ] 

C 0 



In any linear wave refraction model, shoaling and 
refraction coefficients are first computed and then the wave 
height is determined. 

3 . Linear Wave Refraction Model 

Munk and Arthur (1952) described the theory of wave 
intensity along a refracted wave. Dobson (1967) followed 
the theory and developed a linear refraction model. The 
Dobson model is briefly described here. 

The definition sketch for the wave ray equation is 
shown in Figure 2.1 in that A^B]^ is tangential to circular 
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y 

A 




► x 



Figure 2.1 Definition Sketch for Wave Ray 
Equations (Dobson, 1967) 
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arcs of radius C 6 t and (C+5C)6t and S and n are very small. 
Two basic characteristic equations for the propagation of a 
wave ray are the ray curvature equation 



R 



Da 
= DS 



1 DC 
C Dn 



( 2 . 8 ) 



and an equation for the ray separation factor B (Dobson, 
1967) 



dB 2 dB DC 

, . 2 dt DS 

dt 



6 C 




(2.9) 



Dobson manipulated these equations using Gauss 
integral theory into a convenient form where Eqs . (2.8) and 
(2.9) can be easily applied in terms of cartesian 
coordinates and time variable as 



R 



^ [ sin a 



3C 

dX 



cos 




( 2 . 10 ) 




+ P 



dB 

dt 



+ q 6 



0 



( 2 . 11 ) 



where 



*> ri ^ p 

P(t) = 2 [cos a — + sin a 7 — J 

oX oy 



( 2 . 12 ) 



2 

2 " C 

q(t) = C[sin a — ^ - 2 sin a cos 
3x Z 



3 2 C 

3x9y 



2 

2 3 C, 

+ cos a — j] 

° Y 



(2.13) 
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The phase speed is a function of the coordinates so 
that C = C(x,y). The wave ray can be constructed with 

Equation (2.10) within a coordinate system. The ray curva- 
ture R is computed using finite difference methods between 
two consecutive points along a wave ray by iterating 

Act = R A S (2.14) 

at each step. So the new coordinates can be computed as 

X S+'S = X s + AS COS a (2.15) 

y S +AS = y s + AS cos a (2.16) 



where : 



R 



R s + r s+as 



a 



a s + a S+ AS 

2 



The ray path is constructed grid to grid over the sea bottom 
until it reaches the required location. 

The ray separation factor 6 in Equation (2.11) is 
computed applying Fox's formula for a second order linear 
equation. This is a finite difference method and requires 
that the separation factor be known at the starting and 
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preceding points. This condition is conveniently satisfied 
by starting the ray in the deep water where g is equal to 1. 

The characteristics of the wave at any point in the 
coordinate system may be determined from the solution of 
Equations (2.10) and (2.11) provided the speed, C(x,y) is 
known. However it is common to consider waves of different 
deep water characteristics travelling towards the coastline 
over a known bathymetry, so that it is more convenient to 
consider the problem in terms of water depth, d(x,y) rather 
than wave speed. Thus these equations are expressed as 
functions of the mean water depth. 

Dobson used a quadratic least square surface to fit 
depth locally at intermediate points using 

d = Pi + P 2 x + p 3 y + c 4 x 2 + P 5 xy + P 6 y 2 (2.17) 

where the partial derivative of depth, d, with respect to x 
and y is expressed as functions of p. Therefore the local 
water depth can be obtained by interpolation method. 

The bathymetry used in this refraction program was 
obtained from the original NOAA data and projected onto a 
six second (latitude and longitude grid) modified Universal 
Transverse Mercator (UTM) grid after screening the data for 
errors and bad points. Insufficient data were available for 
off-shore bathymetry. Thus it required interpolation to 
fill the empty grid cells. Piece-wise linear triangulation 
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was used to interpolate to a 200 meter rectangular grid. A 
nine-point weighted linear average was applied to smooth 
bathymetry (Oradiwe, 1986; McGee, 1986). 

B. DIRECTIONAL WAVE SPECTRUM 

1 . General 

Directional characteristics of ocean waves are 
important, not only for understanding wave generation 
mechanisms but also for practical problems such as wave 
forecasts in the design of coast and harbor protection 
works. Directionality is equally important in the 
refraction of ocean waves. Various techniques incorporating 
floating buoys, pressure gauge arrays, electromagnetic 
current meters and SAR have been utilized to determine 
directional properties of the ocean waves. In this study 
data have been acquired from the NDBC heave/pitch/roll buoy 
for deep water spectra and pressure gauge arrays for shallow 
water. The theoretical background of these measurements is 
presented in the following sections. 

2 . Deep Water Wave Spectrum 

Longuet-Higgins (1963) suggested a method to compute 
directional wave spectrum from water surface elevation and 
two orthogonal components of wave surface slope. The 
surface elevation is supposed to be the result of the 
superposition of a large number of long-crested progressive 
sinusoidal component waves in random phase (Battjes, 1972) . 
Several other investigators (e.g., Cartwright, 1963? Long, 
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1980; Mitsuyasu, 1975; and Long and Hasselmann, 1979) have 
suggested and refined the technique. 



The directional spectrum S(f,6) as given by Longuet- 
Higgins (1963) is expressed in the form 



where S(f) is the one-dimensional frequency spectrum and 
D(f,a) is angular distribution function with the condition 

2 " 

/ D(f ,a)da = 1 (2 . 19) 

0 

Equation (2.19) can be expanded by a Fourier series 



The coefficients of the first two harmonics can be estimated 
from the co- and quad-spectra from heave (index 1), and 
slope signals (indices (1,2), correspond to x,y slopes) 
(Longuet-Higgins , 1963; Long, 1980) as 



S(f,a) = S ( f ) • D(f,a) 



(2.18) 



oc 




n=l 



( 2 . 20 ) 



2 77 



a n = / D ( f , a ) cos a da = 

0 



K 



( 2 . 21 ) 



2 7T 

bi = / D ( f , a ) sin a da = 

0 



K 



( 2 . 22 ) 



13 



a 2 



(2.23) 



2tt 

= \ D(f,a) cos 2a da 

0 



c - c 

u 22 c 33 



k 2 c 
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b 2 



2 71 

/ D ( f , a ) sin 2a da 

0 



2C 



33 



k 2 c 
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where the constant K is defined as 



(2.24) 



K- [ C22 * C33 ] 1/2 
C 11 

A heave/pitch/roll buoy and pressure sensor array 
provide a measure of n , r\ x , and n y at fixed points. The co- 
and quad-spectra of pairs between these three quantities are 
computed using the Fast Fourier Transform. 

The wave directionalities are determined by the 
directional distribution function. Equation (2.20) may be 
rewritten as 



D(f , a) 



+ — [r, cos (a-a n ) + r 0 cos 2(a-a-)] 
2 TT T, 1 12 2 



(2.25) 



where 



r x = [a 2 + b 2 J 1/2 and r 2 




+ b 



2 1/2 
2 J 



In the case of unimodal spectrum, Longuet-Higgins 
(1963) proposed a cosine powered distribution function 
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(2 . 26) 



D(f,a) 



G ( s ) 



| cos 



1 , 

2 ( a-a ) 



2S 



where 



_ , Q . _ j2S-l r 2 (S+l) = 2S-1 (Si) 2 

G(S) _ r (2S+1) (2S) ! 

is a normalizing function to satisfy Equation (2.19). The 
parameter S is a function of wave frequency and a is the 
mean direction computed from 

. -l. b l, . -1 / Q 13. 

a = tan ( — ) = tan (-=- — ) 

a l g 12 

Cartwright (1963) showed that the value of S is determined 
by 




1 + 3r 2 + (r 2 + 14r 2 + 1) 1/2 
S 2 = 2(1 - r 2 ) 

The cosine distribution model is used most in uni- 
modal wave analysis. The same expression is used in this 
study . 

3 . Predicted Shallow Water Wave Spectrum 

The shallow water wave spectrum can be predicted 
from the deep water spectrum by multiplying the squares of 
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both refraction (K R ) and shoaling (Kg) coefficients and the 
Jacobian of the direction function by each frequency 



component. Hence the deep water spectrum is transformed 
into shallow water spectrum by 



Since the wave period remains invariant, J only depends upon 
the changes of local wave direction. 

4 . Measured Shallow Water Wave Spectrum 

The two-dimensional shallow water wave spectrum is 
determined from the data acquired from shallow water 
sensors. The same technique as described previously for 
deep water spectrum is used. In this study, the pressure 
gauge arrays have been utilized in which subsurface pressure 
fluctuations are measured to infer the surface elevation 
(Kim, 1974; Pawka et al . , 1976). 

From linear wave theory, the relation between wave 
elevation ( r ) and subsurface pressure (P) is described by 



2 2 

S(f,«) = S 0 (f 0 ,a 0 )K R K s J 



(2.27) 



where the Jacobian 




J d (f , a) 



n (t) 



WObU ^ p/4-\ 

cosh k(d + z) J v ; 



cosh kd 



(2.28) 
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where 



d = water depth in meters, 

z = vertical distance from mean sea level to 
pressure gauge (z is positive upward from 
sea level) , 

P = pressure measured in meters of water head 
about a zero-mean pressure. 

The term in the brackets in Equation (2.28) is the spectral 
transfer function relating pressure to surface elevation as 
a function of frequency and depth. Equation (2.28) can be 
modified to include differential distance and elevation 
between two sensors of an array to compute co- and quad- 
spectra of surface slope as 



S, (f) = 

^ AX 2 

dx 



, cosh kd -.2 , f . 

cosh k(d+z) J AP 



S , (f) 
n,dn 

dx 



1 r cosh kd . „ m 

x cosh k(d+z) J T> ,AP x VI ' 



(2 . 29) 



(2.30) 



where Ap x is subsurface pressure difference between two 
sensors on the x-axis and Ax is differential distance along 
the same axis. Similarly, pressure power spectra and co and 
quad-spectra can be computed along y-axis. Once the power 
spectra, co, and quad-spectra are obtained, the directional 
spectral density may be determined using the same technique 
outlined earlier for deep water directional spectrum. 
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III. DATA ACQUISITION AND SENSORS 



A. LOCATION AND TIME 

The area of Monterey Bay was selected as the study area 
which poses complex problems due to the highly irregular 
bathymetry of the Monterey submarine canyon. Three sensors 
were utilized for acquisition of wave data. The locations 
of the sensors are as follows (see Figure 3.1 for bathymetry 
and location) . 



No. 


Tvoe of Sensor 


Location 


Depth (M) 


1 . 


NDBC roll/pitch 
wave buoy 


3 6°48 ' 00 ' 'N 
122°23 ' 59 ' 'W 


2003 


2 . 


Array of pressure 
gauges (4) 

Marina 


3 6°4 2 . 0 'N 
12 1°4 8 . 9 'W 


15.0 


3 . 


Array of pressure 
gauges (4) 

Santa Cruz 


3 6°57 . 0 'N 
122°00.2'W 


13 . 1 



Selected wave cases from January 1988 are used in the 
analysis. During the winter season, weather conditions are 
generally favorable for creating higher energy waves. In 
addition, an abnormal storm was encountered in the region. 
To select a particular date and time, the conditions for 
narrow band wave spectrum were emphasized to simplify the 
interpretations . 
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Figure 3.1 Location of Sensors and Bathymetry in Monterey Bay 



In comparing coincident data sets, on the average a time 
lag of one hour was assumed for the waves to travel from the 
deep water buoy to the shallow water array sites. 

B. DEEP WATER WAVE DATA 

Deep water wave data for the pitch/roll buoy were 
obtained from the NDBC. The NDBC 3-meter buoy measures 
heave, pitch and roll movements with respect to 3-axes 
magnetometer. The magnetometer axes are aligned with fore 
and after (bow positive) , athwartship (starboard positive) 
and vertical (positive upward) axes of a buoy. The analog 
sensor measures vertical acceleration and displacement in 
addition to hull pitch and roll. The magnetometer bow and 
starboard output, together with vertical displacement, pitch 
and roll output from analog sensors are sent to the 
Directional Wave Data Analyzer (DWDA) after voltage 
amplification and filtering. 

The hull azimuth angle is calculated from bow and star- 
board magnetic components with necessary corrections made 
for pitch, roll and hull magnetic fields. The magnetic azi- 
muth is converted into true azimuth using the magnetic vari- 
ation for the site. The true azimuth is then used to 
calculate East-West and North-South components of hull 
slope. One hertz time record of hull displacement and 
slopes are stored in DWDA for spectral analysis and trans- 
mission to shore station through GOES (Steele et al., 1985). 
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Data provided on magnetic tape included environmental 
parameters, wave spectra, directional wave parameters and co 
and quad spectra. Directional wave parameter data were 
utilized to determine directional wave spectrum. The 
unimodality of the directional waves was confirmed by 
studying the regional weather maps to conveniently apply 
cosine angle distribution function. Two-dimensional spectra 
are shown in Figure 3.2. 

C. SHALLOW WATER WAVE DATA 

The calibrated data for arrays of pressure gauges were 
obtained from the Scripps Institute of Oceanography. The 
array consists of four pressure gauges located in a six- 
meter square area which measures pressure in centimeters of 
water. Each sensor is located about 50 to 100 cm above the 
bottom. Pressure power spectrum is calculated which in turn 
is transformed into surface elevation and surface slopes 
n x and ny, spectra by applying linear wave theory transfer 
function as described in Chapter II. 

For the computation of the cross-spectra, the following 
values were utilized: 

- sampling interval = At = 2 sec, 

- total data points = 2048, 

- data points for one sample = 128, 

- total record length = 32 min, 

- Nyquist frequency = 0.25 Hz 
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SPECTRAL DENSITY (M**2-HZ-DEC) 



DIRECTIONAL WAVE SPECTRUM (DEEP WATER) 
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Figure 3.2 Part of Two-dimensional Deep Water Wave 
Spectrum, 3 January 1988, 0800 PST 
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- frequency bandwidth = Af = .00391 Hz, 

- degrees of freedom = v = 32. 

The spectra were cut at frequency of 0.2 Hz in the 
higher frequency range due to spurious overestimation of the 
transfer function. The wave period for high quality 
estimate is approximately 5-20 seconds. 

To ensure accuracy and to check computations, the 
results were compared and verified with the published data 
in the Coastal Data Information Program of January, 1988. 
Two-dimensional spectrum at the peak energy frequency is 
depicted in Figure 3.3 for the Marina array on 29 January 
1988. 

D. OCEAN WAVE TRANSFORMATION PARAMETERS 

The Dobson linear refraction model has been used to 
estimate the shoaling transformation of the wave direction- 
al spectrum. The reverse projection method is employed. 
That is, the refraction program was run backward from 
shallow to deep water, commencing at a shallow water loca- 
tion in required depth, for Marina 15 meters and for Santa 
Cruz 13.1 meters. Rays of a particular frequency were 
propagated offshore in increments of .2 degrees (0.2°) over 
the range of all possible incoming wave angles. Once the 
depth is deep, the rays are stopped and deep water wave 
directions are recorded. The rays were then returned to 
retrace the same path to verify the backward refraction and 
to calculate transfer coefficients Kpj, Kg and J. 
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SPECTRAL DENSITY (M**2-HZ-DEG) 



MEASURED SHALLOW WAVE SPECTRUM-MARINA 




DIRECTION (DEG) 



Figure 3.3 Two-dimensional Shallow Water Wave 
Spectrum Measured at Marina on 29 
January 1988 
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To include refraction over the shelf to the north of 



Santa Cruz, supplemental bathymetric data were meshed into 
the model. The shelf was more coarsely gridded to a 1 Km 
grid because of the large area involved and due to the fact 
that the waves are primarily in deeper water while 
traversing the shelf. The transfer parameters calculated 
from the two areas were multiplied together to determine a 
close relationship for use in spectral transformation. Some 
example ray traces are shown in Figure 3.4. 

1 . Shoaling Coefficient 

. . 2 

Shoaling coefficients, K s , computed from the model 
for Marina and Santa Cruz sites are depicted in Table 3.1. 



TABLE 3.1 

SHOALING COEFFICIENTS, Kg 



Frea. 


Shoaling 


Coefficients 


(Hz) 


Marina 


Santa Cruz 




(15 m) 


(13.1 m) 


0.030 


2 .204 


2.351 


0.040 


1.689 


1.796 


0.050 


1.389 


1.472 


0.060 


1.197 


1.263 


0.070 


1.068 


1.121 


0.080 


0.979 


1.022 


0.090 


0 .918 


0.951 


0.100 


0.877 


0.902 


0.110 


0.851 


0.868 


0.120 


0 .837 


0.847 


0.130 


0.834 


0.836 


0.140 


0.839 


0.834 


0.150 


0 . 850 


0.839 


0.160 


0.867 


0.850 


0.170 


0.887 


0.865 
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FREQ = 0.030 HZ 




X 

FREQ = 0.060 IIZ 




X 

Figure 3.4 Ray Traces from Marina to Deep Water 
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FREQ = 0.09 0 II 7. 




X 

FREQ = 0.130 HZ 




X 

Figure 3.4 (Continued) 
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2 . Refraction Coefficient 



The model computes refraction coefficient, Rr, for 
every 0.2° of nearshore angle. Kr values were averaged for 
every degree of deep water direction. In the linear refrac- 
tion model, an assumption is made that the underwater topog- 
raphy changes are small and gradual but, occasionally, this 
assumption does .not hold in Monterey Bay because of the 
canyon. As a result, the refraction model sometimes gives 
erratic estimates of the Kp values. Therefore, all Kr 
values above three were truncated to the value of three. It 
is also observed that refraction estimates are more stable 
in higher freguencies (short period waves) than lower 
frequencies (long period waves). See Figure 3.5. 

3 . Jacobian (J) 

da 0 

The Jacobian (pp ) is calculated from the change in 
the deep water direction as a function of near shore angles 
for each wave frequency for all possible wave ray direc- 
tions. A plot of olq versus a s (Figure 3.6) shows a slope of 
nearly 45° over most of the angles close to the beach 
normal. The slope tends to increase as the wave ray 
undergoes greater refraction. The Jacobian calculated from 
these values (Figure 3.6) gives similar results, i.e., the 
Jacobian estimates are constant and closer to one when the 
approach angle is within 25°-30° of the beach normal and 
increasingly unstable/erratic at extreme limits. Values at 
the limits may not be useful for computation of directional 
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Figure 3.5 Refraction Coefficients for Waves at Marina 
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Figure 3.6 Plot of Deep Water Approach Directions vs Ashore 
Angles for Waves Arriving at Marina Location 



spectra particularly in the lower frequency range of .07 Hz 
and less. 



The Jacobian appears to behave similar to the 
refraction coefficients. The Jacobian estimates at higher 
frequencies are more stable than at lower frequencies. 
Figure 3.7 shows the Jacobian estimates at different fre- 
quencies for the Marina array location. 

4 . Transformation of Deep Wave Spectrum 

The deep wave spectrum is transformed to shallow 
water wave spectrum using Equation (2.27). Transformation 
parameters calculated from the refraction model are 
multiplied with deep water directional spectrum to obtain 
shallow wave directional spectrum. Figure 3.8 depicts the 
transformation of 0.09 Hz wave for different wave 
directions. The graph of Kg is not shown since at a 
particular frequency it has a constant value for all 
directions . 
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Figure 3.7 Jacobian (J) at Various Frequencies for Marina 
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Figure 3.8 Transformation of Spectrum from Deep Water 
to Shallow Water at Marina for Peak 
Frequency 0.09 Hz 
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IV. RESULTS AND DISCUSSION 



A. REFRACTION MODEL PREDICTIONS 

Transformation parameters, and Jacobian, determined 
by the linear wave refraction model were variable, with a 
tendency to increased variability with increased refraction. 
The lower frequency waves are more refracted than high 
frequency waves. Therefore, these coefficients tend to vary 
smoothly in higher frequency (short wave period) and 
increased in variation with a decrease in frequency. The 
variability of K R and J is greater for the Santa Cruz site 
for which incident waves must undergo significant refrac- 
tion (see Figures 4.1 and 4.2). The characteristics of each 
location and refraction model predictions are described 
here . 

The Marina location is shadowed by Point Pino in the 
Southwest and Point Santa Cruz in the Northwest. High fre- 
quency waves, 0.11 Hz and above, are limited to directions 
from 235°-310° at the array site, whereas low frequency 
waves may arrive from even southerly directions after 
considerable refraction (refer to Figure 3.4). Thus, the 
overall cut-off limits are 235° in the South and 310° in the 
North. The wave rays, at low frequency, are sparse and tend 
to form caustics which result in higher variability for 
estimates of coefficients. The variation of J also depends 
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Figure 4.1 Comparison of Refraction Coefficients; 
Marina and Santa Cruz 
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Figure 4.2 Comparison of Jacobian at 0.13 Hz 
at Marina and Santa Cruz 
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upon the direction of approaching waves in addition to the 
frequency. If the wave ray coincides with, or is closer to 
the beach normal, the J values smoothly vary over the 
directions. As the angle of approach increases, J tends to 
vary with increasingly high estimates. 

The Santa Cruz site poses more complex problems than 
Marina because of its location. The waves approaching from 
the westerly direction makes almost a 90° turn to arrive at 
the shoreline. The low frequency waves travelling north and 
north-eastward may arrive at this location. The ray pattern 
shows (Figure 4.3 a,b) scattered and crossing rays (caus- 
tics) with occasional or a few wave rays from the west or 
higher approach angle at low frequency. This gives very 
large and highly variable estimates of J. The ray pattern 
becomes more regular with increase in frequency. It is also 
noticeable from the ray patterns that the J estimates are 
generally large and more variable compared with the Marina 
values (Figure 4.2). The results of Kp estimates behave 
similarly. 

B. WAVE DATA SELECTION 

Data for wave measurements for the month of January 1988 
were examined. A number of cases were analyzed out of which 
three are described here. The selected days are 4 January 
0800 PST, 18 January at 0200 PST and 29 January at 2000 PST. 
These periods provide narrow band energy density at differ- 
ent frequencies and direction. Figures 4.4, 4.5, 4.6 and 



37 



FREQ = 0.040 HZ 




FREQ = 0.060 HZ 




(a) Low Frequency 



Figure 4.3 Wave Ray Traces at Santa Cruz 
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FPEC = 0.090 H7. 





(b) At High Frequency 



Figure 4.3 (CONTINUED) 
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Figure 4.4 Measured Spectral Density on 18 January 1988, 
0200 PST 
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Figure 4.5 Measured Spectral Density for 29 January 1988, 

2000 PST 
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Figure 4.6 Measured Spectral Density for 3 January 1988, 
0800 PST 



Tables 4.1 and 4.2 provide a comparison of measured 
directional spectra at the three locations for each day. 



TABLE 4.1 

DEEP WATER ENERGY DENSITY WITH MEAN DIRECTIONS 



FREQ 


3 JAN 


0800 PST 


18 JAN 


0200 PST 


29 JAN 


2000 PST 


(Hz) 


DIR 


ENERGY 


DIR 


ENERGY 


DIR 


ENERGY 


( 1/sec) 


(Deg) 


(m 2 -sec) 


(Deg) 


(m 2 -sec) 


(Deg) 


(m 2 -sec) 


0.0300 


173 . 4 


0.000 


283.4 


0.000 


183 . 3 


0.000 


0.0400 


163.8 


0.000 


297 . 3 


0.000 


234 . 3 


0.000 


0.0500 


240.5 


0.000 


265.0 


0.000 


223 . 7 


0.000 


0.0600 


214 . 2 


0.000 


271 . 5 


26.658 


224 . 1 


0.804 


0.0700 


256.6 


1.245 


284 . 6 


72.766 


256.5 


2 . 125 


0.0800 


270.0 


1.569 


292.8 


70.196 


273.5 


6.346 


0.0900 


268 . 5 


1.198 


302 . 2 


38 . 581 


276.6 


8 . 018 


0.1000 


242.8 


1.510 


311.8 


22 . 677 


277.5 


4 .256 


0.1100 


236.7 


2 . 398 


308.7 


12 . 638 


276.6 


3 . 007 


0.1200 


228 . 3 


2 .755 


313 . 5 


11.789 


272 . 2 


2 . 329 


0.1300 


227.7 


2.881 


321.3 


10.475 


268 . 8 


1.570 


0.1400 


223.9 


2 . 190 


313 . 6 


8 . 089 


269.5 


1.113 


0.1500 


217 . 6 


1.699 


321.7 


3 . 894 


268.7 


0.675 


0.1600 


210.1 


1.418 


322.5 


2 . 687 


265.5 


0.609 


0.1700 


190.5 


1.347 


312 . 1 


2.315 


259.3 


0.350 


0.1800 


179 . 2 


0.889 


309.6 


1.940 


259.6 


0.233 


0.1900 


175.8 


0.899 


317 . 2 


1.550 


267 . 2 


0.184 


0.2000 


173 . 6 


0.834 


311.2 


1.147 


245. 1 


0.146 


0.2100 


165.2 


0.670 


313.8 


0.767 


235.6 


0.166 



TABLE 4.2 

WAVE CHARACTERISTICS OF SELECTED CASES 



CHARACTERISTICS 

H sign. 

Freqpeak 

Frec 2mean 
Direct ion mean 



3 JAN 88 
2.2 m 
0.13 Hz 
0.18 Hz 
230° 



18 JAN 88 
7.9 m 
0.07 Hz 
0.10 Hz 
300° 



29 JAN 88 
2.3 m 
0.09 Hz 
0.12 Hz 
275° 
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The California coast was hit by an unusually severe 
storm during 16-18 January 1988 which resulted in high 
waves, strong winds and sea level well in excess of 
predicted values (Cayan, et al., 1988). Widespread damage 
occurred along the Southern California coast. This storm 
developed about 300 NM west of San Francisco on 16 January 
1988 and traveled south-east at about 33 knots. The storm 
center passed off Monterey Bay on the morning of 17 January 
1988 continuing on its SE track and making landfall at Avila 
Beach in Central California. The bay experienced record 
high waves 14 hours after the storm passed from the area. A 
significant wave height of 7.9 meters was recorded on the 
morning of 18 January 1988. The dominant wave period was 
about 14 sec with a direction of 284° with narrow band 
energy measured at deep water and shallow water locations 
( Figure 4.4). 

Deep water energy density with mean direction for 
selected cases are given in Tables 4.1 and 4.2. On 29 
January 1988 measured spectral density of deep water waves 
and shallow coastal site are shown in Figure 4.5. 

Reviewing Tables 4.1 and 4.2 and other published data, 
deep water waves generally travel eastward during the month 
of January with slight variability of about 20° on either 
side of true East. This direction is dominant in the 
energetic middle frequency range of 0.08 to 0.13 Hz. At 
other frequencies, the direction tends to be variable. 
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The third case of 3 January 1988 provides a different 
view from the normal as it contains sea wave energy at 0.13 
Hz frequency with a mean southerly direction of 228° (see 
Figure 4.6). In this case, deep water waves undergo less 
refraction to arrive at the Santa Cruz array. 

C. COMPARISON OF MEASURED AND CALCULATED SPECTRA 

The directional spectra from wave measurement data were 
determined using a cosine-powered distribution function 
(Equation (2.26)). The value of the S parameter, in most 
cases, was less than 20 and relative directional energy was 
in a broad band of 50° to 60° spread. 

Predicted and measured shallow wave directional spectra 
for Marina are compared in Figures 4. 7-4. 9. The deep water 
wave approaching directions are limited between 235° to 
310°. These limits change slightly with decrease in wave 
frequency. In the low frequency range, the predicted 
spectral energy is overestimated (see Figures 4.7 and 4.8). 
The deep water spectra of high frequency waves from the West 
can be transformed into shallow water with reasonably good 
estimation comparable with actual wave measurements. The 
predicted spectrum from westerly waves is shown in Figure 
4.9. The predicted mean shallow water direction tended to 
be more aligned with the deep water spectrum than the 
measured shallow water spectrum. 

The transformation does not appear to appropriately 
shift the mean direction of spectral peak energy of deep 
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Figure 4.7 Comparison of Shallow Water Spectra at Marina 
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Figure 4.8 Comparison of Shallow Wave Spectra at Marina 
for 18 January 1988 
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Figure 4.9 Comparison of Predicted and Observed Shallow Wave Spectra 
at Marina for a Frequency of 0.09 Hz for 29 January 1988 




water wave spectra. Two additional cases of 5 January 1988 
and 30 January 1988 are shown in Figures 4.10 and 4.11 for 
high frequency spectra. The mean approach direction for 
these frequencies in deep water and shallow water are: 

Deep Shallow 



5 


January 


1988 


0.12 


Hz 


255° 


272° 


30 


January 


1988 


0.13 


Hz 


292° 


281° 



The predicted spectrum in shallow water also indicates high 
spectral density in the mean direction. 

At higher frequency sea waves, the refraction model 
overpredicts spectral density by less than 18% when summed 
over all angles at a particular frequency. For swell waves, 
the results are largely exaggerated. 

Results of the Santa Cruz location are depicted in 
Figures 4.12 and 4.13. The results are quite different from 
Marina. The predicted spectra are highly overestimated at 
lower as well as higher frequencies. The relatively simple 
case of 3 January 1988 was particularly selected for Santa 
Cruz to consider ocean waves from the southwest direction. 
Even this simple case does not provide satisfactory results. 

These differences are broadly associated with bathy- 
metry. The non-linear effects due to diffraction and 
currents were neglected which may introduce errors. It is 
surmised that the complexity of bathymetry in Monterey Bay 
may contribute significantly to overestimation and variation 
of transformation parameters. 
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Figure 4.10 Comparison of Predicted and Observed Shallow Wave Spectra 
at Marina at a Frequency of 0.12 Hz, 5 January 1988 
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Figure 4.11 Comparison of Shallow Wave Spectra at Marina 
for 30 January 1988 
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Figure 4.12 Comparison of Shallow Wave Spectra at 
Santa Cruz for 29 January 1988 
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Figure 4.13 Comparison of Shallow Water Wave Spectra at 
Santa Cruz for 3 January 1988 at 0800 PST 



V. CONCLUSION 



Deep water directional wave spectrum was determined from 
data acquired from the NDBC buoy and transformed to shallow 
water using the linear wave refraction model. The 
transformed spectrum was compared with measured shallow wave 
spectra at Marina and Santa Cruz. 

The transfer parameters of refraction coefficients and 
Jacobian predicted by the linear refraction model are highly 
dependent upon wave period and direction of approaching 
waves. The Jacobian appears to be more sensitive to wave 
period and direction than the refraction coefficient. These 
coefficients are generally stable at higher frequencies 
and/or when the wave ray is normal or nearly normal to the 
coastline . 

The model tends to overestimate the transformed 
spectrum. The estimates are less than 18% in excess of the 
measured spectral density for easterly waves at frequencies 
of 0.09 Hz and higher for the Marina site. At lower 
frequencies the results are largely exaggerated. 

The model provides unrealistic large variations in wave 
height calculations for waves at Santa Cruz. This is 
presumably due to its location, orientation of shoreline 
with prevailing north-easterly wave direction and the 
presence of Monterey submarine canyon. 
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It is clear that the strong irregularities of bathymetry 
in Monterey Bay causes non-linear effects. The linear 
refraction model is not suitable for handling such problems. 
Further modeling efforts including non-linear effects are 
needed for Monterey Bay wave study. 
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